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It is a generally known problem that the behaviour predicted from perturbation theory for asymp- 
totically free theories like QCD, i.e. asymptotic scaling, has not been observed in Monte Carlo 
simulations when the series is expressed in terms of the bare coupling go- This discrepancy has 
been explained in the past with the poor convergence properties of the perturbative series in the 
go. An alternative point of view, called Lattice-Distorted Perturbation Theory proposes that lat- 
tice artifacts due to the finiteness of the lattice spacing a cause the disagreement between Monte 
Carlo data and perturbative scaling. Following this alternative scenario, we fit recent quenched 
data from different observables to fitting functions that include these cut-off effects, confirming 
that the lattice data are well reproduced by gQ-PT with the simple addition of terms ff(a n ). 
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1. Asymptotic scaling 

The so-called j3 -function quantifies the dependence of the coupling constant on the lattice 
spacing a: 

Integrating this expression one obtains a relation between g and a which is the usual expression for 
the running of the coupling: 

where 

f PT (g 2 )=e S£F (b g 2 ) 2b ° (l+d 2 g 2 ) (1.3) 

is the 3-loop scaling function, bo, b\ and ^ the usual 1, 2 and 3-loop coefficients, and A is a constant 
of integration. For any lattice prediction of QCD to have physical relevance it should follow the 



asymptotic scaling according to Eqs. (1.2) and (1.3), in the limit of the bare coupling go — > 



Figure 1 is an example of the strong disagreement between Monte Carlo data and go— PT, 
where the data is taken from [[I]]. 1 We fit the same set of data in different renormalized schemes 
known in the literature; the gy schemes were introduced in [0], the g£ in [0], and the gEi in [0]. 
We can see that the mismatch decreases using a renormalized coupling constant instead of the bare 
one [§. 

However, a better reduction (see Fig. 1) is obtained with Lattice-Distorted PT using go as 
expansion parameter and including lattice artifacts ff(a n ) due to the systematic error in the Monte 
Carlo data due to the finiteness of a. The expression (O) becomes: 



-If 2n a l 
a L (So 

where A^ is the scale parameter of lattice QCD 



o) = T - TTT + £ c' n (g 2 )f P " T (g 2 )), (1.4) 
/prfeo) n=l 



2. Lattice-Distorted PT 



Eq. ( |L4| ) can be written in a more useful form by including the leading &(a n ) and next-to- 
leading &(a n ') terms only: : 



fPT{gl) 



G lfpA G l) 



8 v o'fM) 



(2.1) 



,! ' v G V '^(G2) 

where Go is some convenient, reference value of go, e.g. we take Go = 1 for the Wilson data 
(corresponding to j8 = 6). In Eq. (L4), the values of n,n',v and v' 2 depend on the lattice action 



'The lattice spacing in this figure, a rc is introduced in [|l|] and defined via the force, F(r), in an analogous fashion 
to r , but with r^F c (r c ) = 0.65. 
2 Normally v = v'. 
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P 

Figure 1: Plot of data points a -1 obtained from r c with the Wilson action together with fits from next-to- 
leading order lattice-distorted PT and from PT in different renormalised schemes (see text). The data was 
taken from |jlj. 
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Si 
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.13(1) 


-.028(4) 


.27 






K-K* 


4.6(5) 


.22(6) 


-.07(3) 


.60 


DBW2 
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ro 


1.24(3) 


.064(6) 


-.019(4) 


4.1 




i 


T c 


1.248(4) 


.2756(6) 




3.7 



Table 1: Results from the fits of lattice data to lattice distorted PT, using the 3-loop fpj for the Wilson 
action and the 2-loop fpj for the Iwasaki and DBW2 actions . 

being used in the Monte Carlo simulations, and on the quantity being used to set the scale, a (see 
Table 1). Eg. for T c with the Wilson action, we have n = 2, n' = 4, and v = v' = 0. 

The physical quantities used to set a in this study are: the string tension a |Q, [6|]; the length 
scale ro and r c [p], 0|, and the critical temperature T c The K,K* mass point on the plot of My 
versus Mp S , where My and Mp$ are the mass of the vector and pseudoscalar meson respectively 
[^], provides the fifth source of a~ l data [[J. As well as the Wilson action, data from the Iwasaki 
and DBW2 action were also studied. 

Where possible we fit data to next-to-leading-order (i.e. Eq. ([Q|)) obtaining the coefficients 
X n v and Y n iy shown in Table 1. At NNLO the results generally don't change from NLO, so we 
consider the results in Table 1 as our best fits. The only exception to this is the case of K,K* mass 
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Figure 2: Plot of the Monte Carlo data a 1 obtained from observables with the Wilson action together with 
the NLO lattice-distorted PT curves. The 3-loop fpj function was used. 
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Figure 3: Plot of ci d ^ ta /aJ t versus a PT for all data in the Wilson case, where a PT = Al/ fprigo)- Data from 
the NLO fit and the 3-loop fpj function were used. 
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point, where the NNLO fit finds a A/, value much closer to that from other quantities. The DBW2 
fit for T c was constrained to be at LO due to the small number of data points available. The Wilson 
data was fitted to 3-loops, whereas the Iwasaki and DBW2 actions to 2-loops. 

We plot the data and fits for a~ 1 versus j8 in Figures || and [| for the Wilson and Iwasaki actions 
respectively. In order to show the high level of agreement between the data and lattice distorted 
perturbation theory, we plot the ratio of the data to the fit in [3] and || showing agreement at the 
percent level. 

3. Conclusions 

Note that the quality of the fits (as indicated by the % 2 in Table 1, and Figs. 3 & 5) is generally 
excellent. Furthermore, the value for each action is remarkably consistent (see remarks in Sec. || 
regarding the K — K* fit). Note also that the series in a appears convergent, in that Y n iy «X n v 
for all cases. Finally, X n v , which is a measure of the leading-order lattice systematics, is ~ 20% 
in the Wilson case, and, as expected, it is significantly smaller for the improved actions. (Note that 
there are only three T c data points for the DBW2 action, hence the fit parameters in this case can 
be discounted.) 

Converting Al for the r c and T c Wilson case to A^ we obtain 3 

A^ =0 = 217±21 Mev, (3.1) 

MS ' v 

in agreement with previous lattice determinations in quenched QCD [Q, |]]. 

Finally, since the lack of perturbative scaling is probably due to a mixture of both lattice 
artefacts and the poor convergence of the go— PT, we have started to fit data to Lattice-Distorted 

3 We are still investigating the extraction of Aj^ from the Iwasaki and DB W2 fits. 
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Figure 5: Plot of a c { ata /afi t versus a PT for all data in the Iwasaki case, where a PT = Al/ fprigo)- Data from 
the NLO fit and the 2-loop fpj function were used. 



PT using a renormalized coupling constant instead of go (see [§]). In particular, in the g£ scheme 
where the 3-loop coefficient is known, we obtain results consistent with [f|], but the addition of the 
0{a n ) terms improves the quality of the fit. 
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